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We show that the “practical” Bell inequalities, which use intensities as the observed variables, 
commonly used in quantum optics and widely accepted in the community, suffer from an inher¬ 
ent loophole, which severely limits the range of local hidden variable theories of light, which are 
invalidated by their violation. We present alternative inequalities which do not suffer from any (the¬ 
oretical) loophole. The new inequalities use redefined correlation functions, which involve averaged 
products of local rates rather than intensities. Surprisingly, the new inequalities detect entangle¬ 
ment in situations in which the “practical” ones fail. Thus, we have two for the price on one: full 
consistency with Bell’s Theorem, and better device-independent detection of entanglement. 


I. INTRODUCTION 

Since the first derivation of Bell inequality based, ex¬ 
cept for “freedom”, only on the assumption of locality 
and (a form of) realism (or hidden causality, states, vari¬ 
ables) which one can find in the CHSH paper Q, many 
other ones were derived. They either apply to differ¬ 
ent number of settings, or different number of systems 
and observers, or different observables (for a recent re¬ 
view of this “zoo” see Ref. 0)- Here we would like 
to study the case of Bell inequalities which are specially 
tailored for general quantum optical situations, especially 
the cases for fields of varying intensities (essentially, un¬ 
defined photon numbers). The case of strongly driven 
type-II parametric down conversion (with multiple pair 
emissions), which we shall call bright squeezed vacuum 
will be our working example here. But, obviously, mul¬ 
titude of other examples can be given. We use in the 
title the phrase “quantum optical fields” to stress that 
the inequalities we study do not assume a fixed (overall) 
photon number. 

Such inequalities play an important role in singling 
out non-classical phenomena in the optical domain. E.g. 
with the use of such inequalities Bell’s theorem for ap¬ 
proximate EPR states was finally given in 0. The idea 
was to use different observables than the ones discussed 
by EPR . In 0 observables with singular Wigner rep¬ 
resentations were used (essentially, photon number par¬ 
ity operators measured after a displacement). Still, in 
general quantum optics there is an obvious demand for 
device-independent entanglement indicators (witnesses) 
and therefore Bell inequalities for variables straightfor¬ 
wardly related to the usual measured observables. We 
shall propose such. But first we shall present the prin¬ 
cipal work which addressed the problem, and contains a 
kind of mainstream approach. 


A. Standard approach to Bell Inequalities for 
quantum Optics 

We shall review here the derivation by Reid and Walls 
fl], which set the standard in the field, for many years. 
Two observers measure intensities of light at outputs of 


some local interferometric devices. It is customary to 
think of them as outputs of polarization analyzers. We 
shall follow this picture here, as our aim is to convey to 
the reader the general ideas. Let us concentrate on the 
derivation of Clauser-Horne-like Q inequalities, as they 
are often thought of as the most basic ones. 

The approach of Ref. [lj would involve the following 
elements of the hidden variable description. We have in¬ 
tensities potentially measurable by Alice in the two out¬ 
puts of he analyzer Ia ± ( 9 , A), which are non-negative (an 
extension of the approach of to this regime 0 was given 
in Q) and in principle unbounded. For Bob we have 
/ B± (0, A). Here, 9 and <f> symbolize the local settings of 
the polarization analyzers (the can be replaced, if needed, 
by description characteristic for elliptic polarization mea¬ 
surements, like local Bloch vectors, etc.), controllable by 
the observers, and A stands for hidden variables (hid¬ 
den causes not present in the quantum description). The 
hidden variables are assumed to be distributed with a 
probability density p( A). Up to this point all that seems 
to be perfectly in line with the usual thinking behind Bell 
inequalities 

Had the above been the only set of assumptions about 
local hidden variable description, one would not expect 
any loopholes in the resulting inequalities. However, next 
in Ref. [lj comes a remark that the “total intensity 
through each polarized is written” 

Ia(X) = Ja + (M) + Ja_(M), 

Ib (A) = Ib+ (</>, A) + Ib_ (</>, A). (1) 

It is argued that Ia( A) and Ib( A) are in principle inde¬ 
pendent of 9 and <f >, and that they correspond to local 
intensities with polarizers removed. All this is an addi¬ 
tional assumption on the form of hidden variable theo¬ 
ries, which does not have to hold generally in all possible 
hidden variable descriptions. E.g. polarizer in such de¬ 
scriptions could ’’enhance” for some values of A the total 
intensity registered behind it, and for some lower it. Let 
us present an old argument. Note that in classical op¬ 
tics if one has two polarization filters on the way of a 
beam, and the first one transmits the horizontal light 
only and the second one only the vertical one, then no 
light exits such two polarizer-system. However insertion, 
between the two, of a polarization filter selecting diago- 
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nally polarized light results in an output intensity which 
is 1 /4 of the original one. Thus, the presence of polarizers 
does not always reduce intensity. One could here intro¬ 
duce a version of ”no enhancement” assumption, namely 
I A ( A) > Ia + {0, A) + I A _(0,X), etc., but this also is an 
arbitrary constraint on hidden variable models. 

The additional argumentation against such an assump¬ 
tion can come directly from quantum optics. Take as an 
example (bright) squeezed vacuum. In such a state the 
intensity is undefined. However, one should not here that 
the condition perfectly holds for classical stochastic the¬ 
ory of light, and this could have been the intuition that 
led the authors of [I] to adopt it. Still, the condition does 
not hold in the case “stochastic electrodynamics [§]. Even 
if we argue that the condition reflects energy conserva¬ 
tion, this can be challenged by the fact that the field plus 
the system of detectors is in fact an open system. Thus 
the assumption seems not so much justified, afterall. One 
of the aims of this work is to formulate quantum optical 
Bell theorem for systems of varying intensities which is 
free of such assumptions. 

The additional assumption 0 allows one to construct 
local hidden variable correlation functions which are con¬ 
strained by CHSH-like [3j| inequalities. 

m4>) = -(U TT) - <2) 

where 

((i A+ (o) - i A _mi B+ (<i>) - i B -m 

= f p(Wa + {0, A) - I A. (9, X ))(I B+ (</>, A) - I B _ (</>, A)), 

(3) 

while 

(I A I B ) = fdXp(\)I A (\))I B {\). (4) 

with all that one quickly arrives at the usual form of the 
CHSH-like inequality: 

I E(0, </>) + E{0, <fj) + E(9', 0) - E(6', cf>')\ < 2, (5) 

for details see Ref. pj. This possible only due to the 
fact that D = ( I A I B ) is the common denominator in all 
correlation functions entering the CHSH-like inequality. 
Thus assumption 0 is crucial in the derivation. 

In the case of Clauser-Horne-like inequalities one pro¬ 
ceeds similarly. However, in the case of considerations in 
Ref. [I] the no-enhancement assumption is used instead 
of (|T]). Namely that there exists variables, which are 
independent of the local setting (and the remote one...) 
I A (X) and I B ( A) such that 

I A {X) > I A+ (d, X), 

Ib{X) > I B+ {4>, A). (6) 

This allows to introduce functions 

G(M) = J d\p(X)I A+ (0, X)I B+ ((f), A), (7) 


and 

r A {0) = j dXp(X)I B (X)I A+ (6,X), 
r B (4>) = J dXp(X)I A (X)I B+ (9,X). 

( 8 ) 

The variables I A (A) and I B ( A) are treated as intensities 
(potentially) measured by the local observers “without 
the polarizers present”. The additional assumption © 
is not in conflict with 0, as it logically follows from 
it. Still, it wides the class of hidden variable theories 
under consideration. Nevertheless, is prone to the same 
criticism as 0. 

To derive the actual inequality one can use the Clauser- 
Horne lemma: for real numbers 0 < x, x' < X and 0 < 
y, y' < Y one has 

xy + xy' + x'y — x'y' — xY — Xy < 0, (9) 

(in Ref [lj there is a typo in the equivalent inequality). 
Thus 

I A+ (0, X)I B+ ((f, A) + I A+ (0, X)I B+ (</>', A) 

+I A+ ( 6 A )I B+ (0, A) - I A+ {9', X)I B+ (</»', A) 

-Ia + (9,X)I b ( A) - I A (\)I B+ (<i>,\) < 0. (10) 

This after averaging allows one to formulate a Clauser- 
Horne-like inequality 

G(0,^,) + G(6»,^) + G(0 , ,0) 

— G{6', (ft) — r A (0) — r B {4>) < 0. (11) 

The above inequalities, like most of other Bell inequal¬ 
ities, e.g. the original CHSH ones, can be used in the 
analysis of the formalism of quantum optics and its pre¬ 
dictions aimed at proving the impossibility of having a 
local realistic model of all that. Both types of inequalities 
can be violated in quantum optical experiments involv¬ 
ing quantum entanglement. As a matter of fact, in very 
many quantum optical experiments the inequalities were 
used in analysis of results. Still, the validity of such re¬ 
sults (theoretical and experimental) is limited by the fact 
that we use the additional assumptions, already at the 
level of derivation of the inequalities. Other loopholes 
(detection loophole, locality loophole), which are due to 
experimental imperfections, can pop up in the analysis 
of the experimental results, however these are not a trait 
of the above inequalities. 

II. NEW APPROACH 

We shall show that one can derive Bell inequalities for 
quantum optics of a different kind, but still being general¬ 
izations of the CHSH ones and CH ones, which are free of 
the above additional assumptions and thus loophole free 
(the loopholes can now appear only in the methods of ex¬ 
ecution of the actual experiments). We shall additionally 
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argue that they are usually more strongly violated than 
the ones presented above. 

The inequalities will use “rates” instead of intensities. 
We define them as 

^ ± ( 0 ) = W’ ( 12 ) 

= ( 13 ) 

and from the operational point of view they show the 
relative distribution of the measured intensities at the 
two outputs (we still use the two output working exam¬ 
ples, but all that can be generalized to more complicated 
cases) of the local polarizer. We do not assume that the 
total intensity is independent of the local setting. Instead 
we define the following variables 

i b{4>) = E*=± I B i {<t>)- (14) 

(15) 

Of course in quantum theory averages both quantities do 
not depend on on the local settings of the polarizes (well, 
in the experiment one could expect some variation due 
to imperfections affecting losses), but we introduce such 
a description having in mind the hidden variable theories 
for which “anything goes”. 

Let us begin with a derivation of Clauser-Horne-like 
inequalities. 

A. CH-like inequalities 

In this derivation we shall use a hidden variable model 
of the rates. Our only basic tool will be I A± (0, A) and 
/b ± (0, A). In the hidden variable description we shall 
assume existence of values for all possible settings of the 
following variables: 

R ^ =I -m§ <16) 

where I A (d, A) = Ei=± /a<( 0, A), and 

p O n_ J B+ (0. A ) / 17 N 

r b + {& x ) I A ^X) ’ (17) 

with Jg(0, A) = E»=± lBi{<t>, A). Whenever the given de¬ 
nominator vanishes , for certain A and certain local set¬ 
ting, we additionally define that for the given (hidden) 
rate one has Rj ± (a, A) = 0, where J = A or B, and 
a = 6,<j), respectively. 

Obviously one has 0 < Rj + (a, A) < 1. Thus the 
Clauser-Horne lemma can be used to get 

- 1 < Ra + (0 , A )R b+ (0, A) + R a+ (0, A )R b+ (0', A) 

+ R a+ (0', A )R b+ (0, A) - R a+ ( O' , A )R b+ (0', A) 

-R A+ (e,X)-R B+ {^X)<0. (18) 


This after averaging over p( A) leads to 

- 1 < K{9, 0) + K{9, 0') + K(6', 0) - K(9', 00 

-5 a+ (0,A)-5 b+ (0,A) <0, (19) 

where 

K(0, 0) = J dXp(X)R A+ (0, A )R b+ (0, A), (20) 

and 

S J+ (a) = J dXp(X)R J+ (a,X). (21) 

This is the new form of a loophole free Clauser-Horne 
type inequality for quantum optical experiments. The 
loopholes are moved to the experimental imperfections, 
like inefficient detection, etc. 

The reader might be worried that we divide here by the 
intensity, and how should look like the quantum defini¬ 
tion for equivalent expression. Of course this depends on 
the model of detection used, which defines the intensity 
operators at the exists of the local devices (polarizers), 
to illustrate this we shall use the simplest model: the 
photon number operator. 

As we have four outputs, and an example involving po¬ 
larizations, let us introduce annihilation (and creation) 
operators a B , ay, b B and by (creation operators would 
have a dagger, f). Here H and V stand for horizontal 
and vertical polarizations, and a, b denote propagation 


beams, 

towards Alice and Bob, respectively. The 

opera- 

tors for arbitrarily oriented (linear) polarization analyz- 

ers are 

cos 0a B + sin day = ae. 

(22) 


— sin 9a B + cos 9ay = ag ±, 

(23) 


cos 4>b B + sin (f>by = bp , 

(24) 


— sin 06jj + cos (j>by = bp±. 

(25) 

Let us 

assume the the operators with _L denote — 

exit of 

the local polarizers, and the other ones + output, 
the chosen intensity model given by 

Under 


C^Cq;, 

(26) 


lJ - , 

(27) 


where c = a, 6, a = 6, 0 for J = A, B, respectively, the 
quantum versions of the functions in the inequality read: 

K(e,ct>) Q =Tr[eR A+ (6)R B+ ((l>)\, (28) 


and 

S J+ {a) Q =Tr[gR J+ (a)}, (29) 

where g is the quantum state, and Rj + (a) are given by: 

R J+ (a)=ilj° Ij+ flj°. 

l^l = ± *J± 


( 30 ) 
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The operators IIj 0 are projectors into the subspace, of 
the Fock space, of the inodes described by annihilation 
operators ch and cy, which does not contain vacuum. 

In case the notation seems a bit too dense, let us write 
an example of such an operator explicitly. E.g., the op¬ 
erator R a+ (0) acts in the joint Fock space of modes ag 
and a g ±, and reads 

Ra + { 0 ) = ni° t Q *°; -ni° (3i) 

a g a g + a g± ag ± 

with 

ni° = 4 - |0,0) aa (0,0|, (32) 

where I a is the identity operator in the space, and 
|0,0) a is the vacuum of the space satisfying ae|0,0) a = 

a 9 r|0,0) a = 0. 

As in quantum optics violations of Bell inequalities of 
the type presented in jlj usually decreases with the grow¬ 
ing number of photons, see e.g. [16], the inequalities in¬ 
troduced here might be less prone to this effect, because 
their form warrants that that terms with higher photon 
numbers contribute less to correlation functions A(0, 0)q 
and to the local averages Sj + (ct)Q. Thus we shall ben¬ 
efit from their use twice. They are not only formulated 
without additional assumptions (which is a kind of built 
in loophole), but also we shall gain on the sensitivity of 
the inequalities in detecting non-classicality. This will 
be shown in the examples of the next Section. But be¬ 
fore this we shall derive CHSH-like Bell inequalities for 
quantum optics based on similar concepts. 


while this is a correct inequality, its usefulness for quan¬ 
tum optics is highly limited. This is because in most Bell 
experiments the produced states are such that the vac¬ 
uum components dominate. Especially this is the case for 
parametric down conversion, see e.g. [1]. As a result the 
values of the correlation functions are very low, and the 
left hand side cannot breach the bound of 2. Still such 
inequalities could be of use in the case of event ready 
experiments, see 0- 

To devise CHSH inequalities which are immune to this 
deficiency, one can proceed as follows. One can introduce 
new rate functions which behave differently when the in¬ 
tensity is zero. Simply one can redefine the (arbitrary) 
value of Ra + (0, A) and Rb + (0, A) for zero total intensity, 
and put is as 1. No change is required for Rj_ ( a , A). We 
can denote the redefined rate functions by R' A ( 6 , A) and 
R' b (0, A). Again we have 

—1 < Rj + (a, A) — Rj_ (a, A) < 1. 

With the new variables the hidden variable model for the 
correlation function reads: 

<0(0,0) = J d\p(X)[R' A+ (9, A) — R a _ (9, A)] 

x [R'b+ (0; A) — Rb- (0, A)]. (37) 

This definition does not suffer from the mentioned defi¬ 
ciency, because in the case of vacuum the values 1 pop 
up. We have a modified inequality: 

I C(9, 0) + <7(0,0') + <7(0', 0) - <7(0', 0')l < 2- (38) 


B. CHSH-like inequalities 


The CHSH-like inequalities can be derived using 
the same basic hidden variable objects I A± (0, A) and 
0B±(0, A), and their functions R A+ (9, A) and R B+ {4>, A), 
plus additional ones 


i?A_ (0, A) = 


7a(0,A) 


Ia{9, A) 

where I A {0, A) = X)»=± I Ai (0 , A), and 

/a(0, A) 


Rb- (0, A) = 


04(0, A) ■ 


(33) 


(34) 


Recall the all R variables equal zero for zero total lo¬ 
cal intensity. Obviously this can be done by introducing 
correlation functions 

F(9, 0) = J dXp(X)[R A+ (0, A) - R A _ (0, A)] 

x [Rb+ (0, A) — R b _ (0, A)]. (35) 

As — 1 < Rj + (a, A) — Rj_ (a, A) < 1, one can proceed in 
the usual way to derive: 


|F(0, 0) + F(9, 0') + F(0', 0) - F(0', 001 < 2. (36) 


In the case of vacuum state the modified inequality, fol¬ 
lowing from the above definitions, is saturated , while for 
some other quantum state it can be violated, see further 
down. 

We still have to give a quantum formula for C(0, 0)q. 
To this end one must define the operator R'j (a). This 
will be again done for the trivial model of intensity in 
terms of the number operator. 

R'a + (0) = Ra + (0) + |0, 0) aa (0,0| (39) 

and 

^B + (0)=^B + (0) + |O,O) 6fc (O,O|, (40) 

where 10, 0)& is the vacuum for the b modes. 


III. EXAMPLES 

Surprisingly the new inequalities, not only are free of 
an inbuilt loophole (that is the additional condition), but 
one can easily find examples of quantum optical states 
with violate them more robustly than the old-type in¬ 
equalities. We shall consider here our working example 
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of the four mode (bright) squeezed vacuum. The state 
reads 


I BSV) = -j- Y, \/^+ltanh"r|^ n) ) I (41) 
cosh r 


where 


ivl n) > 

= Em=o(- 1 ) m l n - m >«n \m)a v \m)b H | n - m ) bv , 

(42) 

where a and b refer to the two directions along which 
the photon pairs are emitted, H/V denote horizon¬ 
tal/vertical polarization, and T represents an amplifica¬ 
tion gain. The state is rotationally invariant, i.e., all 
probabilities are dependent only on 0 — ip, and this can 
also holds for elliptic polarizations.. 

First, let us compute the value of the CHSH expression 
in fashion of Ref. [I|, according to formula IQ. For 
every Fock component of state (HTl) with 2 n photons , 
the correlation function is a cosine function with one of 
its minima at 9 — <p = 0 and amplitudes |(2 n + n 2 ). 
Thus, the numerator in Eq. m, summed over n > 0 

with weights —“> reads 

{{i A+ {e) - i a _(6)){Ib + {cP) - i B _m) 

= — 2 cosh 2 F sinh 2 T cos(0 — p), (43) 

while the denominator is equal to 

( IaIb) = — 2cosh2r + ^cosh4T + cosh 4 r. (44) 


If we optimal settings 9 = n, O' = |7r, p = and <p' = 
— j7r , we get violations of the Bell inequality for T < 
0.4911. 

In the new loophole-free approach utilizing the ratios, 
for an individual 2n-photon component we have 

(^- ) I [R'a + W - Ra _ (0)] [R'b + (P, A) - R b _ (P, (n)> 

= -^eos^- 0), (45) 

while <O,O,O,O|[^ + (0) - R A -miR' B+ (<t>^) - 

Rb_ (p, A)]|0, 0,0, 0) = —1. After averaging over Fock 
components with weights (n + 1) tanh 2rl (r) cosh -4 F, we 
get 


C(0,P)q = [ x - I cos (9-p) (4 log cosh r 

+(3 + cosh 2 F) sinh 2 r)] . (46) 

This put into the CHSH expression provides violation 
for T < 0.8866. Not only does inequality m provide a 
fully valid Bell inequality for quantum optical fields, it 
also extends a region of gain T, for which they can be 
violated. 



FIG. 1. Interference contrast (visibility) for quantum pre¬ 
diction for K(6,P)q (new) and G(6,P)q (old). The usual 
threshold for violation of CHSH ineqialties is given by the 
straight line, just to indicate the scale of improvement. 


A. CH-type inequality 


After some algebraic manipulations the coincidence 
rate function K(9,P)q (l28l) for the bright squeezed vac¬ 
uum state (TUI) can be found to be 

m P) Q = 4(1- ^Ar) + (47) 

x (—13 + 12 cosh2r + cosh4r + 8 In cosh T). 

The visibility (contrast) of the interference is therefore 


Vnew (T) 


Kmax(.@ •> 4*) Q Kmini®i 0 ) Q 
Kmax(@i 0 )Q T Kmini ®? 0 )Q 
3 + cosh 2 T - 4 lnc ffip r 

_ I _Sinh- r 

3 + 3 cosh 2 T 


(48) 

(49) 


where the minimum and maximum is taken over p and 
9. In the case of the appraoch of [lj the visibility is given 
by (see e.g. 0) 


v 0 id( r) 


i 

1 + 2 tanh 2 T 


(50) 


The visibilities for rates are much more pronounced, see 

Fig. ffl 


IV. EXTENSION OF THE METHOD 

We have also compared the thresholds for these two ap¬ 
proaches when adopted to the chained inequalities with 
L settings per side. Such inequalities (for the new ap¬ 
proach) read 

\C(6 1 ,p 1 ) + <7(02,00 + C(6 2 , p 2 ) 

+ ...+C(Q L ,p L )-C{e u p L )\<2(L-\). (51) 

This was done for equal spacings of the settings 0i = 
^ 2 * 4 L ^ an d c t )i = t£- Reid-Walls-like chained inequal¬ 
ities are violated for values of gain T decreasing with 
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L 

yycrit 

1 R-W 

y\crit 
-L ratios 

L 

j^crit 

1 R-W 

j^crit 

J- ratios 

2 

0.491 

1.229 

8 

0.251 

1.425 

3 

0.408 

1.342 

9 

0.237 

1.431 

4 

0.355 

1.372 

10 

0.224 

1.437 

5 

0.318 

1.392 

11 

0.214 

1.441 

6 

0.290 

1.406 

12 

0.205 

1.445 

7 

0.268 

1.416 

13 

0.197 

1.448 


TABLE I. Critical values of Fr_w and F ra tios above which 
the Reid-Walls-like and ratios-based chained Bell inequalities 
cannot be violated. 


which we still observe the violation of the inequalities are 
presented in Table HU As expected, chained inequalities 
for ratios are violated for higher gains T than the origi¬ 
nal ones. A simple explanation is that in the latter, the 
the contributions from high numbers of produced pairs 
can grow unboundedly, while in the ratios-based version, 
they cannot be higher than 2. Notice, for example, that 
to violate the inequality for ratios at T = 1.4 one needs 
only 4 local settings, whereas it takes 7 settings to violate 
the original inequality. 

V. CLOSING REMARKS 


L 

-perit 

Chained 

-perit 

Chained, ratios 

2 

0.915 

1.123 

3 

1.053 

1.367 

4 

1.165 

1.482 

5 

1.260 

1.586 

6 

1.345 

1.687 

7 

1.427 

1.795 


TABLE II. Critical values of gain rgh ained and Teamed, ratios 
above which chained inequalities for moduli of intensity dif¬ 
ferences and and similar ones for ratios cannot violated in 
function of the number of local settings L 


growing L. Surprisingly the behavior of the ratios-based 
inequalities is opposite!. Table H] lists the critical values 
of Tr-w and r rat i os , below which the chained inequalities 
of “length” L are violated. 

The new method gives improvement also in the case of 
chained inequalities recently derived in Ref. [TlJ . The in¬ 
equalities, if one uses the wording of this paper, are based 
on averages of \Ia + (9) — Ib + (<(>)|, which can be shown to 
have properties of geometric distance. Let us replace this 
by (\Ra + (9) — i?B + (</>)|). The resulting inequalities read 

ELi(\RA + m-R B+ tti)\) 

+ (\Ra + (9i+i) ~ R B+ {4>i) |) 

> {\Ra + (0i)-Rb + {cI> l )\). (52) 

We have numerically studied the violation of the 
chained inequalities with the cutoff of at most 25 pairs 
produced in each run. The values of the gain, below 


One can generalize the approach to situations involving 
more observation stations, compare Ref [k|, and Bell 
experiments involving generalizations of Bell inequalities 
for many parties, like 171. 

The new approach in fact leads to a possibility of a new 
look at Stokes parameters for non-classical light. This 
will be presented in a forthcoming paper E3 , where it 
will be shown that with the new approach we can get 
better entanglement witnesses (conditions). This is done 
by replacing the Stokes parameters in old entanglement 
conditions of the type presented in Ref. (Till [hit by 


c' <a\ _ R_LO a e ae a S_ L ae ± A_L0 
RAW) — H A —f-T-• 


agag 


a 9 I a ®i 


(53) 


As a matter of fact all entanglement witnesses for states 
of two-qubit systems can be directly mapped into entan¬ 
glement indicators, for optical four mode fields (of the 
kind studied here), involving polarization measurements. 
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